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ABSTRACT 

A  revised  model  for  predicting  surface  bubble  loss  (SBL)  as  a  function  of  wind 
speed,  surface  grazing  angle,  and  frequency  is  discussed.  SBL  is  defined  as  the  resid¬ 
ual  loss  (i.e.,  the  loss  in  excess  of  that  due  to  spreading  and  chemical  absorption)  in 
acoustic  energy  propagating  near  the  sea  surface  that  is  attributable  to  attenuation 
by  near- surface  bubbles.  This  loss  is  a  key  input  parameter  for  predicting  bubble  at¬ 
tenuation  in  simulations  of  weapon-system  performance  in  near-surface  environments. 
The  original  SBL  model  was  introduced  by  the  Applied  Physics  Laboratory,  Univer¬ 
sity  of  Washington,  in  1993.  The  revision  consists  of  a  new  wind-speed  threshold  for 
breaking  waves  and  the  subsequent  onset  of  bubble  production,  and  an  alternative 
formulation  capable  of  handhng  more  general  ray  geometries  such  as  a  ray  vertexing 
near  the  surface.  Model/data  comparisons  and  examples  of  the  model  appfied  to  a 
linear  sound-speed  profile  are  presented. 
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EXECUTIVE  SUMMARY 

In  the  report  APL-UW  TR  9307,  a  model  was  introduced  for  the  loss  of  acoustic 
energy  propagating  near  the  sea  surface  that  is  attributable  to  bubbles.  This  loss 
is  referred  to  as  the  surface  the  bubble  loss,  or  SBL.  This  report  documents  (1)  a 
revision  in  the  wind-speed  threshold  given  in  the  original  model  for  breaking  waves 
and  subsequent  onset  of  bubble  production,  and  (2)  an  alternative  formulation  of  the 
SBL  model  capable  of  handling  more  general  ray  geometries  such  as  a  ray  vertexing 
near  the  surface.  The  alternative  formulation  consists  of  a  depth-dependent  bubble 
attenuation  coefficient. 

The  revised  model  is 

SBL(dB)  =  [/i.57^o.85^  U>6m/s  (1) 

smd 

=  SBL1„^,  U  <  6  m/s  , 

where  U  is  wind  speed  measured  10  m  above  the  sea  surface,  /  is  acoustic  frequency 
in  kilohertz,  and  6  is  the  nominal  grazing  angle  of  the  surface  bounce  path  {0  >  0°). 
In  the  revised  model,  the  wind-speed  threshold  has  been  increased  from  4  m/s  to 
6  m/s  and  the  form  of  the  decay  for  wind  speeds  <  6  m/s  has  been  changed  to  reduce 
the  likelihood  of  predicting  anomalously  high  SBL  at  very  low  wind  speeds. 

The  alternative  formulation  of  the  model  is 

/*10  (y(z') 

SBL  (dB)  =  2  /  .  ,  (2) 

Jz„  sm[^(z)J 

where  a{z)  is  the  depth-dependent  bubble-attenuation  coefficient,  Zy  is  the  turning 
point  depth  of  the  ray,  and  6(z)  is  the  local  grazing  angle  of  the  ray  at  depth  2r. 
The  depth-dependent  bubble-attenuation  coefficient  is  integrated  between  a  depth  of 
10  m  and  Zy  to  form  an  estimate  of  the  two-way  depth-integrated  bubble  attenuation, 
which  is  equivalent  to  SBL.  An  equation  is  given  for  a{z)  as  a  function  of  wind  speed 
and  frequency. 

Starting  with  Eq.  (2),  the  report  derives  useful  analytic  expressions  for  SBL  based 
on  a  hnear  sound-speed  profile.  Unless  the  ray  undergoes  vertexing  near  the  surface, 
Eq.  (1)  gives  the  same  result  as  Eq.  (2),  and  the  user  should  consider  this  fact  in 
deciding  whether  to  employ  the  simpler  version  in  Eq.  (1)  vs  the  depth-dependent 
version  in  Eq.  (2).  For  ray  vertexing,  Eq.  (1)  is  no  longer  valid,  and  the  formulation 
given  in  Eq.  (2)  must  be  employed. 
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1.  INTRODUCTION 

In  the  report  APL-UW  TR9307,^  APL-UW  introduced  a  model  for  the  residual 
loss  (i.e.,  the  loss  in  excess  of  spreading  and  chemical  absorption)  of  acoustic  energy- 
propagating  near  the  sea  surface  that  is  attributable  to  bubbles.  Since  the  bubbles 
are  generated  by  breaking  -waves  and  reside  within  a  few  meters  of  the  sea  surface, 
this  loss  is  referred  to  as  the  surface  bubble  loss,  or  SBL.  This  report  documents 
(1)  a  revision  in  the  wind-speed  threshold  given  in  the  original  model  for  breaking 
waves  and  the  subsequent  onset  of  the  production  of  bubbles,  and  (2)  an  alternative 
formulation  of  the  model  capable  of  handling  more  general  ray  geometries  such  as 
a  ray  vertexing  near  the  surface.  The  alternative  formulation  consists  of  a  depth- 
dependent  bubble- attenuation  coefficient. 

Section  2  presents  the  revised  wind-speed  threshold  and  includes  some  model/ data 
comparisons.  Section  3  presents  the  alternative  formulation  of  the  model  and  includes 
examples  of  its  use  with  a  hnear  sound-speed  profile.  A  summary  and  notes  on  model 
usage  are  presented  in,  Section  4. 


2.  MODEL  REVISION 

The  revised  model  is 

SBL(dB)  =  1:26^x^  ^1.57^0.85^  U>6mls  (1) 

sin  y 

=  SBLlu^,  [/  <  6  m/s  , 

where  U  is  the  wind  speed  measured  10  m  above  the  sea  surface,  /  is  the  acoustic 
frequency  in  kilohertz,  and  6  is  the  nominal  grazing  angle  of  the  surface-bounce  path 
{6  >  0°).  The  model  includes  a  wind-speed  threshold  of  6  m/s  for  the  production  of 
breaking  waves  (also  known  as  the  Beaufort  velocity)  and  the  subsequent  production 
of  bubbles.  The  revised  model  shown  in  Eq.  (1)  differs  from  the  original  model  only 
in  the  wind-speed  threshold  and  in  the  functional  form  of  the  decay  for  wind  speeds 
less  than  the  threshold  speed.  The  wind-speed  threshold  has  been  increased  from 
4  m/s  to  6  m/s  to  reduce  the  likelihood  of  predicting  anomalously  high  SBL  at  very 
low  wind  speeds. 
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Figures  1-3  show  representative  model/data  comparisons  for  data  collected  be¬ 
tween  20  and  40  kHz  in  the  open  ocean  from  the  research  platform  FLIP^  and  in 
coastal  waters  off  Whidbey  Island.^  Because  the  data  were  collected  at  different  graz¬ 
ing  angles,  all  data  have  been  normalized  to  a  10°  grazing  angle  using  the  formula 

SBLaO”)  = -^  SBLW  , 

sm  10° 

where  SBL(^)  is  the  actual  measured  loss  at  surface  grazing  angle  6.  The  model  (Eq. 
(1)  with  $  set  to  10°)  provides  a  reasonable  representation  of  both  the  open-ocean 
and  coastal  data,  with  nearly  all  the  data  points  falling  within  ±3  dB  of  the  predicted 
curves. 


MODEL  VS  DATA,  20  kHz 


Fig.  1.  Comparison  of  model  predictions  and  measured  20-kHz  data  from  the  FLIP  and 
Whidbey  Island  experiments.  Data  from  various  grazing  angles  have  been  normalized  to 
a  10°  grazing  angle.  Upper  and  lower  dotted  lines  represent  a  ±3  dB  interval  about  the 
predicted  curve. 
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MODEL  VS  DATA,  30  kHz 


Fig.  2.  Comparison  of  model  predictions  and  measured  30-kHz  data  from  the  FLIP  and 
Whidbey  Island  experiments.  Data  from  various  grazing  angles  have  been  normalized  to 
a  10°  grazing  angle.  Upper  and  lower  dotted  lines  represent  a  ±3  dB  interval  about  the 
predicted  curve. 
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MODEL  VS  DATA,  40  kHz 


Fig.  3.  Comparison  of  model  predictions  and  measured  40-kHz  data  from  the  FLIP  and 
Whidbey  Island  experiments.  Data  from  various  grazing  angles  have  been  normalized  to 
a  10°  grazing  angle.  Upper  and  lower  dotted  lines  represent  a  ±3  dB  interval  about  the 
predicted  curve. 
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3.  ALTERNATIVE  FORMULATION  OF  THE  SBL  MODEL 

In  this  section,  an  alternative  formulation  of  the  SBL  model  is  presented  which 
is  capable  of  handhng  more  general  ray  geometries,  including  rays  vertexing  near  the 
surface. 

3.1  Depth-Dependent  Formulation 

The  formulation  in  Eq.  (1)  represents  the  two-way  depth-integrated  bubble  at¬ 
tenuation  for  a  surface-reflecting  ray  under  isovelocity  conditions.  It  is  simple  to 
implement,  and  summarizes  in  a  compact  way  the  effects  of  wind  speed,  frequency, 
and  surface  grazing  angle  on  SBL.  Since  isovelocity  is  assumed,  however,  the  model 
cannot  account  for  refraction  effects,  and  in  the  case  of  ray  vertexing  (9  — >  0)  it 
becomes  invalid. 

By  assuming  an  equivalent,  depth-dependent  bubble- attenuation  coefficient  that 
decays  exponentially  with  depth,  we  can  write  Eq.  (1)  as 

SBL(dB)=2/‘°  .‘"Lf -.-.dz  ,  (2) 

Jzv  sm[0(z)j 

where  a{z)  is  the  depth-dependent  bubble-attenuation  coefficient, 

a{z)  =  ,  (3) 

Zy  is  the  turning  point  depth  of  the  ray,  and  9{z)  is  the  local  grazing  angle  at  depth 
z.  The  ray  reaches  a  turning  point  by  either  reflection  or  refraction;  in  any  case,  z,, 
may  fall  within  the  range  0  <  <  10  m.  We  assume  SBL  is  negligible  for  rays  that 

undergo  vertexing  at  depths  greater  than  10  m,  or  Zy  >  10  m. 

In  Eq.  (3)  Lb  is  the  decay  constant,  or  e-folding  depth,  for  the  exponential  decay 
in  bubble  attenuation  with  depth  and 

ao  =  0.63  X  10-^f/'"7°®VLB,  U>6m/s  (4) 

=  «oUe  t/<6m/s. 

Several  works^~®  incorporate  an  exponential  decay  model  to  describe  bubble  con¬ 
centration  vs  depth  which  can  be  directly  related  to  Ci{z).  A  simple  model  describing 
Lb  vs  wind  speed  for  bubbles  within  the  resonant  size  range  most  relevant  to  torpedo 
frequencies  is  Lb  =  0.07  x  U.  This  equation,  which  describes  measurements  reported 
in  Ref.  7,  produces  estimates  of  Lb  that  fall  approximately  between  estimates  derived 
from  similar  models  reported  in  Refs.  3  and  4. 
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The  attenuation  coefficient  a{z)  is  integrated  between  a  depth  of  10  m  and  Zy  and 
multiplied  by  2  to  estimate  the  two-way  depth-integrated  bubble  attenuation,  which 
is  equivalent  to  SBL.  For  example,  if  the  ray  reflects  from  the  sea  surface,  Zy  equals 
0;  if  the  ray  vertexes,  then  Zy  is  the  vertex  depth.  For  practical  implementation,  the 
maximum  vertex  depth  is  set  to  10  m  because  bubble  attenuation  will  be  extremely 
weak  beyond  this  depth.  An  expression  similar  to  Eq.  (2)  is  given  in  Ref.  8,  which 
discusses  a  model  for  bubble  loss  at  frequencies  of  1-10  kHz.  The  user  is  cautioned 
that  Eq.  (2)  has  a  singular  point  when  6{z)  — >  0  which  will  cause  numerical  problems 
if  not  handled  correctly. 

3.2  Examples  Using  a  Linear  Sound  Speed  Profile 

For  an  arbitrary  sound-speed  profile,  Eq.  (2)  must  be  implemented  numerically; 
however,  instructive  analytic  results  are  possible  using  a  hnear  sound-speed  profile  of 
the  form 

c{z)  =  CQ  +  gz  ,  (5) 

where  cq  is  the  sound  speed  at  the  surface,  g  is  the  sound-speed  gradient  in  inverse 
seconds,  and  ^  is  depth  measured  from  the  sea  surface.  For  a  hnear  sound-speed 
profile,  all  ray  paths  can  be  described  by  circular  arcs  with  a  radius  of  curvature  Rc 
equal  to  (|^|cr)“^,  where  the  ray  launch  parameter  a  is  equal  to  cos[^(2)]/c(2),  with  0 
being  the  local  grazing  angle  at  depth  2  (Fig.  4).  Because  a  is  an  invariant  property 
of  the  ray,  it  can  be  evaluated  anywhere  along  the  ray  path.  For  example,  if  the  ray 
is  surface  reflecting,  a  convenient  point  to  evaluate  a  is  at  the  surface  such  that  a 
equals  cos(0s)/co,  where  Os  is  6  evaluated  at  the  surface.  If  the  ray  vertexes  {6  0) 

near  the  surface,  then  cr  equals  l/c„,  where  Cy  is  the  sound  speed  at  the  vertex  depth 

Zy. 

The  different  cases  of  >  0  and  S'  <  0  are  now  analyzed.  Note  that  s  =  0  means 
the  sound  speed  is  isovelocity  in  the  upper  10  m,  the  ray  paths  are  straight  fines,  and 
there  is  no  advantage  in  using  Eq.  (2)  instead  of  Eq.  (1).  We  begin  with  the  case  of 
downward  refraction,  or  s'  <  0?  ^-nd  the  instance  of  a  surface-reflecting  ray  (Fig.  4A). 
Using  the  small-angle  approximation  for  the  sin  of  the  angle  (valid  for  9  ~  20°),  we 
can  write  Eq.  (2)  as 

rlO 

SBL  =  2ao  /  -^dz.  (6) 

Jo  0[z) 
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Fig.  4.  Ray  geometries  for  a  linear  sound-speed  profile  where  9  is  the  local  grazing  angle, 
9s  —  9  sX  the  surface,  and  ^lo  =  0  at  10-m  depth.  Geometries  in  (A)  and  (B)  show  downward 
refraction  {g  <  0),  with  (A)  showing  a  surface-reflecting  ray  and  (B)  showing  a  ray  vertexing 
at  depth  Zy.  Geometry  in  (C)  shows  upward  refraction  {g  >  0)  and  a  surface-reflecting  ray. 


Now,  again  using  the  small- angle  approximation,  we  can  relate  the  depth  z  to  local 
grazing  angle  9{z)  by 

-  =  f  Wf  -  €)  .  (7) 

recalling  that  6{z)  will  always  be  >  0,  for  <  0.  We  recast  Eq.  (6)  using  9  as  the 
integration  variable,  giving 

SBL  =  2aoRce'^^"  ,  (8) 

J0S 

where  k  =  RcI{2Lb)-  The  limits  of  integration  are  now  from  Os  to  ^lo,  where  Ow  is  0 
evaluated  at  10-m  depth.  Evaluating  this  integral  gives 

SBL  =  ,  (9) 
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where  #  is  the  tabulated  probability  integral,  or  error  function,®  defined  as 


Although  the  solution  is  exact,  evaluation  of  the  difference  between  error  functions 
that  both  have  large  arguments  is  prone  to  numerical  error.  An  asymptotic  evaluation 
of  this  difference,  which  is  more  reliable,  gives 

|$(vTe,„)  -  $(vT«.)] « -^=L,  (10) 

and  thus  Eq.  (9)  reduces  to  simply 

SBL(dB)  =  ^^^.  (11) 

Os 

If  we  replace  Og  with  sin(0),  where  6  is  the  ray’s  grazing  angle  at  the  sea  surface,  this 
is  precisely  the  formulation  for  SBL  given  in  Eq.  (1).  Thus  the  effect  of  downward 
refraction  on  a  surface-reflecting  ray  is  accounted  for  by  simply  using  the  surface 
grazing  angle  9  —  dg  in  Eq.  (1). 

Next  we  examine  a  ray  vertexing  within  10  m  of  the  sea  surface  for  the  same  case 
of  <  0  (Fig.  4B).  In  this  instance,  we  have 

flO 

SBL  (dB)  =  2ao  /  ^7-^  dz  .  (12) 

JZy 

Following  the  above  procedure,  we  get 

SBL  =  2aoRce-"^^^^  (13) 

Jo 

The  solution  of  the  integral  now  gives 

SBL  =  2aoi?ce-"''/^^^y^$(v^dio).  (14) 

The  argument  y/kOw  will  in  many  cases  be  >  2,  making  $  nearly  unity.  Thus  it  is 
sufficient  to  write  _ 

SBL  (dB)  =  2aoy/RcLBTrl2  (15) 

as  the  surface  bubble  loss  for  the  vertexing  ray. 
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For  a  vertexing  ray,  Eq.  (2)  does  give  a  different  solution  than  Eq.  (1),  as  the 
latter  is  invalid  when  0  — >■  0.  Figure  5  shows  SBL  vs  wind  speed  for  a  ray  that 
vertexed  2  m  below  the  sea  surface  owing  to  a  sound-speed  gradient  of  g  =  —0.2  s“^. 
This  value  represents  the  maximum  negative  gradient  within  the  upper  10  m  of  ocean, 
as  measured  during  the  FLIP  experiment^  in  late  afternoon  CTD  casts.  Note  that 
Eq.  (15)  provides  a  useful  guide  for  the  effect  of  Zy  on  SBL.  The  main  influence  is  the 
factor  {Rc  changes  very  little  with  z„).  If,  for  example,  Zy  changes  from  2  to 

3  m,  SBL  is  reduced  by  a  factor  of 

Finally,  the  case  of  upward  refraction,  or  >  0,  and  a  surface-reflecting  ray 
(Fig.  4C)  produces  the  same  result  as  a  surface-reflecting  ray  for  g  <  0  (Fig.  4A) 
although  the  intermediate  steps  are  different.  Ray  vertexing  within  10  m  of  the 
sea  surface  is  not  possible  when  <  0.  These  results  depend  only  on  the  grazing 
angle  being  sufficiently  small  that  the  small-angle  approximation  holds  and  not  on 
the  magnitude  of  g.  For  large  angles,  the  additional  path  length  incurred  (compared 
with  straight-line  propagation)  owing  to  refraction  within  10  m  of  the  sea  surface  is 
small,  and  therefore  Eq.  (1)  remains  a  useful  approximation. 


Fig.  5.  SBL  vs  wind  speed  (frequency  =  30  kHz)  for  a  ray  undergoing  vertexing  2  m 
below  the  sea  surface,  g  =  —0.2  s“^. 
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4.  SUMMARY  AND  ADDITIONAL  NOTES  ON  MODEL  USAGE 

Two  changes  were  presented  in  the  model  for  surface  bubble  loss  (SBL)  introduced 
in  APL-UW  TR9307;  (1)  a  revision  in  the  wind-speed  threshold  for  the  production  of 
bubbles,  and  (2)  an  alternative  formulation  based  on  a  depth- dependent  attenuation 
coefficient  due  to  bubbles,  a{z).  The  change  in  the  wind-speed  threshold  affects  SBL 
predictions  for  wind  speeds  less  than  6  m/s.  The  alternative  formulation  of  the  model 
is  capable  of  handhng  total  refraction,  or  vertexing,  of  a  ray  within  10  m  of  the  sea 
surface,  and  also  any  refraction  that  surface-bounce  paths  may  undergo. 

Useful  analytic  expressions  for  SBL  based  on  the  depth-dependent  formulation  of 
Eq.  (2)  were  derived  using  a  linear  sound-speed  profile.  For  surface  reflection,  it  was 
shown  that  the  model  in  Eq.  (1)  gave  the  same  result  as  the  formulation  in  Eq.  (2). 
The  user  should  consider  this  fact  in  deciding  whether  to  employ  the  simpler  version 
in  Eq.  (1)  vs  the  depth- dependent  version  in  Eq.  (2). 

For  ray  vertexing  (i.e.,  6  0),  Eq.  (1)  is  no  longer  vahd,  and  simple  expressions 

for  SBL  as  a  function  of  the  ray’s  vertexing  depth,  and  the  ray’s  radius  of 
curvature,  were  derived  on  the  basis  of  Eq.  (4).  Because  of  the  exponential  decay 
of  bubble  attenuation  with  depth,  SBL  is  set  equal  to  0  dB  for  rays  that  undergo 
vertexing  at  depths  greater  than  10  m. 

Some  remarks  on  model  usage  follow: 

1.  The  coefficients  in  the  model  are  based  on  acoustic  measurements  made  at  20- 
50  kHz.  The  model  is  therefore  a  high-frequency  model;  as  a  broad  guideline,  it 
is  apphcable  at  frequencies  from  10-100  kHz.  However,  its  results  have  recently 
been  shown  to  be  consistent  with  those  of  a  model  developed  at  the  Naval 
Undersea  Warfare  Center  for  lower  ship-sonar  frequencies.^” 

2.  Scattering  from  bubbles  will  prevent  SBL  from  becoming  arbitrarily  large.  A 
nominal  limit  for  SBL  is  30  dB,  based  on  scattering  from  the  underside  of  a 
uniform  layer  of  near-surface  bubbles.  However,  such  a  limiting  value  for  SBL 
has  not  been  observed  in  data,  because  the  patchiness  of  the  bubble  layer  wiU 
nearly  always  allow  some  surface-scattered  energy  to  pass.  Thus  we  recommend 
an  upper  limit  of  SBL^ax  equal  to  15  dB.  This  value  is  based  on  the  SBL 
measurements  reported  in  Refs.  1,  2,  and  11. 

3.  The  effect  of  near-surface  bubbles  on  the  real  part  of  the  sound  speed  (known 
as  the  sound-speed  anomaly)  is  ignored,  since  within  the  model’s  applicable 
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frequency  range  bubbles  alFect  mainly  the  imaginary  part  of  sound  speed,  or 
attenuation.^^ 

4.  The  model  is  based  on  measurements  from  a  single  surface  bounce.  To  compute 
SBL  for  propagation  channels  with  multiple  surface  (and  bottom)  interactions, 
the  total  SBL  is  considered  to  be  the  sum  of  the  individual  SBL  values  computed 
for  each  surface  interaction.  We  caution  that  the  model’s  accuracy  for  these 
conditions  has  not  been  evaluated. 

5.  The  model  does  not  address  geographical  dependencies.  We  expect  that,  in  gen¬ 
eral,  littoral  waters  and  more  pristine  oceanic  waters  will  differ  in  terms  of  bub¬ 
ble  properties.  With  few  exceptions,  however,  bubble  measurements  reported  in 
the  literature  show  within-region  variability  that  is  greater  than  between-region 
variabihty  (compare,  for  example.  Refs.  2  and  7).  This  fact,  together  with  the 
data  shown  in  Figs.  1-3  for  both  coastal  and  oceanic  environments,  supports 
use  of  a  single  model  at  the  present  time. 

Additional  remarks  concerning  basic  assumptions  employed  in  the  model  are  con¬ 
tained  in  APL-UW  TR9307.^ 
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